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We solve analytically the Langevin dynamics of the classic spherical model considering the ferromagnetic 
exchange and a long-range antiferromagnetic interaction. Our results in the asymptotic regime, shows an equiv- 
alence in the functionality of the spatial and self-correlations between this model and the recently studied 
Ginzburg-Landau frustrated model within the Hartree approximation. A careful discussion is done about the 
low temperature behavior in the context of glassy dynamics. The appearance of interesting features regarding 
the establishment of the ferromagnetic phase is also analyzed in view of the effects of the spherical restriction. 
We propose a new variant of the spherical model in which the global restriction is substituted by an infinite set 
of restrictions over finite size regions. This modification leads to a new dynamical equation that suggests the 
appearance of the low temperature phase transition even in the non-frustrated case were the classic spherical 
model fails. 
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I. INTRODUCTION 

The study of the fundamental role of microscopic in- 
teractions on the relevant properties of magnetic materi- 
als involves two main approaches: discrete models, like 
the Ising model which is commonly used to represent sys- 
tems with uniaxial anisotropy, and continuous models among 
which the Ginzburg-Landau model is one of the most widely 
known. These two approaches have been extensively used 
both numerically! 1-4] and analytically 15.-;8J to describe equi- 
librium and dynamical features observed in real systems. ||9|- 
[12I] Concerning to the local magnetization bound, unlike the 
Ising model in which the site magnetization take values ±1, 
the Ginzburg-Landau double well potential is the responsi- 
ble to guaranty the non-divergence of the local magnetization. 
Unfortunately, this potential introduce a strong non-linearity 
in the dynamical equations describing the time evolution of 
the local order parameter, and do not allow the magnetiza- 
tion of the system to saturate at high enough external fields. 
On the other hand, for discrete models some dynamical treat- 
ments are forbidden due to the natural impossibility of define 
temporal derivatives of the local magnetization. 

In this context, the spherical model ifisi [Till is something 
in between. It is a linear model that can be treated analyti- 
cally, while the saturation of the magnetization is guarantied 
by means of the spherical constraint. Despite this niceties, this 
model fails in obtaining of the ferromagnetic transition in two 
dimensions. On the basis of this disagreement it has been ar- 
gued that the spherical constrain allows the system to be in a 
large number of disordered configurations. In two dimensions 
this degeneracy prevents the system to reach a feiTomagnetic 
phase transition at any non-zero temperature. To manage this 
difficulty a lot of work have been carried out making some 
radical modifications as, for example, the substitution of the 
short range exchange interactions by a long range ferromag- 
netic one. Ill [Tell 

Recently, a great effort has been done in the study of the 
system with competing isotropic interactions, i. e. systems 



whose fluctuation spectrum has a minimum in a non-zero ra- 
dius spherical shell in the inverse space. lfl7l - l2lll In magnetic 
systems, these interaction are usually represented by the fer- 
romagnetic exchange interaction and the anti-feiTomagnetic 
long-ranged dipolar interaction. The Langevin dynamic 
of a Ginzburg-Landau Hamiltonian with such interactions 
have been studied previously by means of the Hartree self- 
consistent field approximation. (a |2ll] In these works the au- 
thors confirm the theoretical picture found analytically by 
means of numerical simulations. Moreover, the symmetry 
breaking properties of the model were also studied by means 
of a very similar procedure. ||22ll In particular, it was shown 
that taking the infinite time limit it was possible to encounter 
the phase transition to the modulated ordered state by means 
of purely dynamical calculations. 

In the present work we solve the Langevin dynamics of 
the spherical model considering the common square field- 
gradient attractive term plus a repulsive long-range antiferro- 
magnetic term in the long times limit. We show that, starting 
from a disordered initial condition, quenches to the T — 
regime develops a long time dynamics very close to that en- 
countered in Ref. |21]. As it is expected, for this frustrated 
spherical model, the value of the critical field dividing the two 
different dynamical regimes is also the saturation value. For 
T > we found a standard paramagnetic relaxation after cer- 
tain characteristic time. Anyway, the divergence of this time 
for low temperatures as well as the existence of modulated 
domains of typical length can be related to the reminiscence 
of a glassy region not completely verified for our model. The 
monotony of the h-T curve for a fixed magnetization is also 
an interesting finding that connects the results with recent ex- 
perimental works. Finally, we propose a new variant of the 
spherical model in which the global restriction is substituted 
by an infinite set of restrictions over finite spatial regions of 
the system. This lead to a new dynamical equation that is also 
presented in detail. We discuss the way in which this modi- 
fication suggest the appearance of the low temperature phase 
transition even in the non-frustrated case were the spherical 
model fails. 
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The rest of the paper is organized as follows. The next sec- 
tion is completely devoted to the presentation of the formalism 
and the solution of the dynamical equation for all the regimes 
involving the value of the applied external field. Some com- 
monly studied observables as spatial and self-correlations are 
presented and discussed in some detail in section |III] compar- 
ing its asymptotic expressions with those recently obtained for 
the Ginzburg-Landau model. In section |IV] the more relevant 
findings are discussed in the light of the known phenomenol- 
ogy of related models and some experimental observations, in 
order to better understand the advantages of using the spher- 
ical model in the context of competitive interactions. Also 
a new dynamical equation is obtained by means of suitable 
modifications to the standard model and is discussed the way 
in which this new approach may lead to the appearance of the 
phase transition. We summarize the most important results 
and discussions in section |V] In appendix |A] the spatial cor- 
relations of some particular configurations are calculated by 
means of geometrical considerations. 



II. MODEL AND FORMALISM 

We consider the following Hamiltonian dependent on the 
local parameter (p{x, t) 

1 



{\I(j)Y ~ 2h(j}{x) 



-—.II d'^xd'^x'J{x-x')(j){x)cj){x')(l) 
2o Jo. Jn 

where J{x — x') represents a repulsive, isotropic, competing 
interaction and Q, is the d-dimensional volume of the system, 
which is finally taken as infinity. The homogeneous external 
field is represented by h and the parameter 5 measures the 
relative intensity between the attractive and repulsive interac- 
tions. The local field is subjected to the restriction 



02(f) d'^x = w?Vt 



(2) 



that represents a system with local magnetization m. In order 
to include this restriction in the dynamics it is widely used the 
effective Hamiltonianll23ll 



[0(f)] =H [0(f)] 



\{t) 



(t>'^{x)d'^x-m^n 



(3) 



and in consequence the Langevin equation of motion is given 
by 

90(f. t) _T , 1 



dt 



+ h- \{t)(t>{x) + -q{x,t) 



) - - / dVj(f-f')0(f') 



(4) 



the solution of this equation is carried out in the inverse space, 
defined by the following Fourier transforms 



(5) 



0(fc) = / d'^x e-**^-"0(f) 



The equation of motion in the momenta space is given by 



d4>{k,t) 
dt 



1 



J{k) + \{t) 



4>{k,t) 



+ h{k)+r]{k,t) 



(6) 



and its solution can be written in terms of the response func- 
tion R{k, t, t') as 



4>{k,t) = (f>{k,0)R{k,t,0) + / R{k,t,t')r]{k,t')dt' 

Jo 

R{k,t,t')h{k)dt' (7) 

/o 



where 



/^ _ ^(^') „-l(fe)(t-t') 



R{k,t,t') = ^e 



(8) 



It has been also defined A{k) = A{k) - A(fco), Y{t) = 
elo i{t')dt' ^ j(^) ^ ^(fc^) ^ ^(fc) = fc2 + 1 

At this point the problem has been formally solved, except 
for the fact that the parameter X{t) remains unknown. It can 
be shown that the spherical restriction leads to the normal- 
ization condition C{t, t) = m^, where C(t, t') represents the 
self-correlation function, obtained from the usual correlation 
function given by 



Cik,k',t,t') = {<j){k,t)(t){k',t')) 



(9) 



Thus, to impose the normalization condition for the self- 
correlation function allows to obtain the following equation 



m^Kit) = A/(t) + 2T / K{T)f{t - T)dT 



where K{t) ^Y{t)2 and 
fit) 



ifl/2(^)g-A(0)(t-r)^^ 



d'^k e-^^e^)* 



(10) 



(11) 



In systems with isotropic competing interactions the fluctu- 
ation spectrum A{k) has a minimum in a shell in the in- 
verse space of certain radius fco. As it was shown in previ- 
ous works[8, 21] the long time behavior of f{t) is dominated 
by the local topology of the fluctuation spectrum around fco. 
In consequence, it is enough to consider an expansion of the 
fluctuation spectrum A{k) up to second order around its min- 
imum, such expansion has the form 



A{k)^^{k- 



(12) 



where A2 = ^^-t^ . With these ingredients it is straight for- 
ward to obtain the following asymptotic result 



(13) 
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where C is certain uninteresting constant depending on sys- 
tem dimension. At this point we carry out the solution of the 
problem separately, first we analyze the zero temperature case 
and then the non-zero temperature case. 



A. Case T = 

In order to solve equation ( flOl i at zero temperature, we 
found suitable to define the function 

S{t)= rA'i/2(r)e-'4(o)(*-r)^^ (14) 

"'0 

and then to split up the original equation into the following set 
of equations 



dSjt) 
dt 



-A{0)S{t)+K{ty^ 



(15) 



with the initial conditions K{0) = 1 and S{0) = 0. Note that 
this system can be put in the form 



^ = -AiO)S{t) + 1 {Am 



(16) 



So, we have to find the asymptotic solution of the above dif- 
ferential equation. In order to do this we compare the terms 
f{t) and S'^{t). In general we can see that there are several 
possibilities varying h. For small fields f{t) dominates the dy- 
namics and in consequence S'^{t) oc /(t). On the other hand, 
fields higher than certain critical value he leads to a different 
behavior. In what follows we obtain the different solutions for 
S{t) and K{t) varying the external field h. 



1. h < he 



2. h — he 

We start by writing equation ( fTSI l in the form 
""•"--^(OlSW + iSari + ^l (22) 



dt 



h^S^it) 



In the long time limit Af{t) ^ h'^S'^{t), then it is possible to 
write 

considering an expansion up to first order of the non-linear 
contribution. Note that if we take now h = he — TnA{0), the 
first order contribution of the non-linear term cancels and we 
have 



dS{t) _ A fit) 



dt m2i(0) 25(i) 
whose solution yields 

A 



/(T)dT 



and 



'TTt ^ I 



(24) 



(25) 



(26) 



Expressions (125 b and (1261 ) has no lower integration limit. This 
is related with the fact that, in the long time limit, any constant 
can be neglected in comparison with the increasing functions 
S{t) and K{t). 



In the case of small enough fields we can obtain that 
S^{t) = A/(t). This consideration implies the relation 

- {Af{t) + h^S^t)) ^ = - (A + h^\) ^ (17) 
m TO 

This leads to the simplified equation 

^ = -A{Q)S{t) + 1 (A + h^Xy f{t)'^ (18) 
dt TO 

whose asymptotic solution is given by 
toA(O) 

At this point to consider the proposed ansatz implies that 

A 



A 



w?A{QY - h^ 



and 



K{t) = XAiOffit) 



(20) 



(21) 



It is worth to note that from the definition of A one can see the 
existence of a critical field he = mA{0). 



3. h> he 

In this case, from equation (fT6l l in the long time regime, it 
is fairly easy to see that 



S{t) oc e(-^(°)+^)* 



and in consequence 



Kit) oc e2(-^(o)+^)* 



(27) 



(28) 



B. Case T > 

At finite temperature the equation ( fTOl i contains a convolu- 
tion term with the function f{t), which leads always to a ex- 
ponential long time behavior, as have been obtain previously 
in literature. In order to clarify this point, we propose the 
ansatzlCjjl 



S{t)^ K^/\T)e'^^°^'-'-^^dT ^^K^t) (29) 
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for the asymptotic dynamic. By doing this, equation ( fTOl i turns 
in an integral equation whose solution can be found by means 
of the Laplace transform. Using expression ( fT3b for f{t) we 
get 



(30) 



where b and ^ are the solutions of the following system of 
equations 



2Tcr(i; 



&2 + 2A{Q) 



(31) 



Finally, it is worth to note that the above system of equations 
always give rise to positive values for both, b and ^. This 
means that at any finite temperature the model reaches the 
same asymptotic dynamics. 




FIG. 1: Dependence of the mean value of the local field {(j)) with the 
external field h for several temperatures. 



III. OBSERVABLES 

At this point we have been able to find the asymptotic be- 
havior of the function K{t) for any value of temperature and 
external field. This allows the calculation of some commonly 
studied statistical observables. In particular, we obtain in this 
section the mean value of the local field (p (magnetization), 
the self-correlation function and the spatial correlation func- 
tion. For some of these magnitudes closed expressions are 
found in terms of the function K{t). 



As can be seen from figure [T] at zero temperature and 
h < he the system behaves like a purely paramagnetic ma- 
terial with a completely linear relation between magnetization 
and external field. On the other hand, for any finite tempera- 
ture the relation between the mean value and the external field 
is no longer linear in such a way that {(j)) is always an increas- 
ing function of the applied external field. 



B. Self-correlation function 



A. Mean value 

Starting from expression (jTJl we obtain the following ex- 
pression for the mean value 



where 



t-)-oo K{t)2 



(32) 



(33) 



This expression implies that at T = 

(0> = 



(34) 



TO for h > he 



while for any finite temperature £^{h, T) is given by the solu- 
tion of the system (ISTl i. In order to clarify such a dependency 
of the mean value with the external field and temperature we 
define the following non-dimensional magnitudes 



he ' A(0)5to2 ' TO 

which are the variables uses in the figure[T] 



(35) 



In order to look inside the ordering process we calcu- 
late now the self-correlation function, defined as C(t,t') = 
{4>{r, t)4i{r, t')) - {4>{f, t)) (0(f, t')). In terms of the function 
K{t) this magnitude can be written 



C{t,t') = 



= / 



^K{t)K{t' 
2T 

^K{t)K{t') 7o 



t + t' 



' K{T)f ( -T]dT 



In this way at T = we obtain 



' (4ft')4 

(*+*')' 



Cit,t') oc < 



(i+t')5 (4ti')3 



for h < he 
for h = hr 



(36) 



-(A(0)+A^)(t + t') 

T for h > hr 



Considering our definition of C {t, t'),it can be seen that for 
high enough fields (h > he) the system presents an interrupted 
aging and eventually evolves to a ferromagnetic configuration 
in which there are no fluctuations at all. On the other hand, 
for h < he the system shows a slow relaxation typical of a 
coarsening scenario. 
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For T > the long time behavior [li^t' > 1) of the self- 
correlations is given by 



t+t' 



(37) 



showing an interrupted aging and an asymptotic temporal 
translational invariant relaxation C{t,t') = C{t — t') typi- 
cal of systems in thermal equilibrium. This means that for 
any finite temperature the system always reach the equilib- 
rium state. 



C. Spatial correlation function 

With the aim of study the texture of the system under dif- 
ferent conditions, we calculate the spatial correlation function. 
In our case the isotropy of interactions and initial conditions 
ensures an homogeneous mean value of the order parame- 
ter. The spatial correlation function is defined as C{R,t) = 
{(j){x,t)(j){x + R,t)) ~ {(j){x,t)){(j3{x + R,t)). At zero tem- 
perature (T ~ 0) we obtain 



C{R,t) 



A 



K{t) J {2ttY 
and assuming koR^ 1 it is possible to obtain 
fit) A 



(38) 



C{R,t) cx 



K{t) (koR)^ 

/ (d-l)7r' 
cos koR - I e 



(39) 



where ko and A2 are defined by the expression ( fT2b . At zero 
temperature the relation f{t)/K{t) contains the dependence 
with the external field. For h < he this functionality tends to 
a constant value, while in the other cases we have 



fit) 

m 



-A{0) + 



Tk)* 



— r- 

t2 



for h 



for h > he 



(40) 



It can be noticed that for fields equal or higher than the crit- 
ical one the system evolves to the completely ferromagnetic 
configuration. This can be appreciated in the decay of the 
spatial correlation function. On the other hand, for h < he the 
presence of the exponential term reveals a phenomenology of 
growing domains with typical length L{t) oc t^^ very com- 
mon in models with non-conserved order parameter These 
domains represent random oriented modulated structures with 
frequency ko. It is worth to note that the power law decay do 
not imply the presence of structural defects necessarily, it is 
possible to show how such a term arise from the average over 
the ensemble of possible configurations (see Appendix lAb. 



When T > it is possible to obtain in the equilibrium state 

d'^k 



C(R) = lim 2T 1 ——re 



ik-R 



^d^k 



{2n) 



-2A{k){t-T) 



„ik-R 



(41) 



where the last integral have been calculated before in Ref. ^ 
and gives 

C{R) cx 3__e-f cos {koR - V') (42) 

(fcoi?) 2 



with 



(d-l)7r d-\ 



tan 



1 



(43) 



Expression (l42l i describes a long times state in which a mo- 
saic of modulated structures establishes with a wave length 
271 /ko and a typical dimension equal to the correlation length. 
As can be seen from ( |43] | this correlation length diverges as 
temperature tends to zero for h < h^ and decreases when the 
field is augmented. It is worth to note that for h> he the spa- 
tial correlation as has been defined tends to zero when T 
indicating an evolution to the homogeneous state. 



IV. PUTTING RESULTS IN CONTEXT 

The results exposed above can also be viewed in the frame 
of the general scheme currently established in literature. From 
a static point of view, the infinite time configuration of the 
system shows no long range stripes order at any finite temper- 
ature. This is consistent with the isotropy of our initial condi- 
tions and the fact that there is no symmetry breaking. In fact, 
what we do is to study the long times dynamic of the ensemble 
of disordered quenches. Anyway, as can be seen in appendix 
IaI the lack of long range order in this case is not due to the 
ensemble degeneration but to the appearance of randomly ori- 
ented domains of typical size. 

On the other hand, it has been reported in literature 1 20l l24ll 
that below certain temperature the Ginzburg-Landau frus- 
trated system behaves as a typical glass. Even when the spatial 
correlation found here at finite temperatures (|42] | is consis- 
tent with the glass configuration, the most basic glassy prop- 
erty, that is, the aging behavior of the self-correlations, is not 
verified in the long times limit. Instead, as we have already 
pointed, there exists interrupted aging for all non-zero tem- 
perature, so the system develops an aging dynamics up to a 
finite time Tg that diverges only at T = 0. This characteris- 
tic time can be obtained by fitting the numerical integration of 
equation ( fTOb with the asymptotic exponential behavior and is 
shown as a function of the temperature in figure|2l We can see 
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0.35 



FIG. 2: Characteristic time Te up to which the aging dynamics can be 
observed. Upper curve corresponds to a field h = 0.5hc and lower 
curve corresponds to h = he- The inset represents the logarithmic 
plot of the main figure. 



From numerical data of figure[T]it has been constructed fig- 
ure [3] in such a way that the plotted curves divides the h-T 
space in two regions. In the upper(lower) region magnetiza- 
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from the inset that the characteristic time behaves as a power 
law in a vicinity of T = 0. As it have been suggested,|25] this 
behavior remains very close to that of the glass in the sense 
that one can always find a temperature below which the aging 
dynamic establishes up to a time far beyond any laboratory 
scale. 

It is worth to note that the divergence of the correlation 
length ^ of the ordered domains at T = have been reported 
in a number of glassy systems. |24i] The spatial correlation 
(3% certainly accounts for such a glassy configuration. As 
it is demonstrated in AppendixlAl this correlation functional- 
ity corresponds to an ensemble of perfect striped configura- 
tions every one of which has long range order in a particular 
direction. This T = spatial correlations do not show any 
more the typical exponential decay but a potential decay term 
that should not be straightforwardly interpreted 
as a direct evidence of the presence of topological defects for 
d > 1. Correlations obtained by means of the average over 
the ensemble of perfect modulated structures (i.e. considering 
all the possible orientations) must be in general a decreasing 
function of R for d> 2, even when we have long range order 
in a particular direction (see AppendixlAl). For the T = case 
the glassy behavior is reinforced from the fact that equilibrium 
is never reached in any finite time. 

Yet, there is another arising feature of this model making it 
interesting for future studies in the context of magnetic frus- 
trated systems. If we go back to figure [1] it may be noticed 
the non-trivial fact that lower temperatures requires lower 
field values to reach the same magnetization. In our spheri- 
cal model approach the completely ferromagnetic phase can 
be observed by means of a total saturation {(f)) = 1 reached 
only at T=0. These features may be suggesting that, for sys- 
tems described with this approach, the ferromagnetic configu- 
ration must be reached at low temperatures by applying lower 
fields. This interesting result has been recently observed in ex- 
perimental systems [121, |26t] and is opposed to the commonly 
accepted phenomenology for the Ginzburg-Landau model. 



FIG. 3: External field as a function of temperature for a given value 
of the magnetization. Dotted curve in the main figure corresponds to 
{(f}) = 0.9 while the inset corresponds to ((ji) = 0.1. The continuous 
lines are fits to quadratic functionalities. 

tion is greater(lower) than certain value (0) > 00 ((</>) < ^o) 
with 4>o = 0.9 for the main figure and ipo = 0.1 for the in- 
set. The monotony and functionality of curves in figure [3] is 
independent of the specific value of (f)o even in the vicinity 
of (/)o = 1 that represents the limit of the modulated struc- 
tures. In other words, if we define as ferromagnetic a more 
real not-completely-saturated state where (0) > 1 — e, then 
figure [3] would be in disagreement with theoretical results for 
other thin film models. |5] Since e can be a sufficiently small 
value, the last definition is pretty near to what can be observed 
experimentally. 

This phenomenology is in good agreement not only with 
experimental works[26] but also with some analytical ap- 
proaches of ultrathin magnetic films considering only striped 
and uniform phases. ll29ll However, as can be seen from the fig- 
ure, we obtain a quadratic fit for the temperature functionality 
that is valid for all values of instead of the more complex 
scaling functionality experimentally observed|12| that is de- 
pendent on the critical temperature. 

So far, it is not easy to understand where this phenomenol- 
ogy encountered in the frustrated spherical model begins to 
be disturbed from the fact that this model does not recover the 
phase transition, even in the pure ferromagnetic case in d = 2. 
Nevertheless, there is an acceptable agreement with the stan- 
dard features established in literature for frustrated systems 
and, moreover, the fact that the local field is bounded in some 
way seems to have good implications in modeling the behav- 
ior of the system in the presence of an external field at T > 0. 
This suggests the necessity of finding a near model in which 
the fundamental limitations already mentioned could be over- 
come. 
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A. Extension of the model 

Any modification of the standard spherical model must be 
done in a way that improves the basic characteristics involved 
in the probable establishment of the ordered phase. Lets con- 
sider again the spherical restriction (|2]i but now writing it in 
the local way 



(44) 



where n is a discrete d-dimensional counter vector and T^j rep- 
resents certain d-dimensional volumes of finite constant size. 
These finite volumes do not overlaps each other and com- 
pletely fill the spatial region il, so the limit i7 — > oo implies 
an infinite number of such finite regions. It is over these finite 
volumes that the spherical constrain is now performed. This is 
a closer approximation to what actually happens in many real 
system. 

Following the same Langevin formalism as in section (HJi 
but now including the multiple restrictions ( l44b we reach to 
the dynamical equation 



t) 

dt 



- T / d'^x'J[x-x')<i>{x') 
Jn 



h - \fi{t)4){x) + ij{x,t) 



(45) 



that is virtually similar to that of (|4|l with the exception that A 
is now also dependent on the space region through the counter 
fi. It means that there is an infinite set of Lagrange multipliers 
each one corresponding to a single finite region V^j. 

It is completely plausible to expect that the particular ge- 
ometry of these finite volumes do not implies substantial dif- 
ference in the phenomenology of the model. So we may take 
them as d-dimensional cubes with a linear size 2a. In this as- 
sumption it is possible to write the equation of motion in the 
inverse space as 



d(l){k, t) 

dt 



1 



+ 7 J(fc) 



(l){k,t) + h{k) + r]{k,t) 



+ E^«W / -0^Hk',t)Qn{k^k') (46) 



where 



)n{k-k') = / e'^^-^"^-^d'^x 

d 



n 



2sin[(/c^ - ka)a] 
k'^ ka 



(47) 



and f— (2n + l)a. 

The more complex but yet linear model thus modified is 
a very good candidate for observing the transition to the or- 
dered phase in lower dimensions. The modification is sub- 
stantial with respect to the standard spherical model in which 
temperature can always destroy the saturated phase, driving 
the system to multiple non-homogeneous configurations and 



yet holding the spherical constraint. In this new approach in- 
finite constraints guaranties that, below certain temperature, 
once the correlation length becomes greater than the typical 
size a of the spherically constrained volumes, these volumes 
behaves as fully-magnetized entities not allowing internal in- 
homogeneities. This fact may eventually drive the system to 
the Ising universality class or, at least, to a new interesting 
phenomenology in which the statistical behavior must have 
a change around some temperature value dependent on the 
physically meaningful parameter a. 

It is worth to note that equation ( |46] | admits the trivial para- 
magnetic solution as the limit case in which all Lagrange mul- 
tipliers are identical Ap = A. As can be demonstrated 



k') = {2TiY5(k^k') 



(48) 



and then equation ( |46] | becomes the already solved equation 
( fTOl ). This limit is thus associated in the modified model with 
the high temperatures region. In this conditions the correla- 
tion length is so small that the spherical restriction over the 
complete system has the same effect that if it were done over 
the finite volumes. 



V. CONCLUSIONS 

The Langevin dynamics of the spherical model has been 
solved in the long times limit considering isotropic compet- 
ing interactions. The general asymptotic functionality of spa- 
tial and self-correlations is found to be equivalent to that of 
Ginzburg-Landau model in the Hartree self-consistent field 
approximation.|21| The work was done studying the ensem- 
ble of disordered realizations, and we do not verify any stan- 
dard glassy transition as that predicted for near models. ll20ll 
Nevertheless, for low temperatures the system evolves to a 
mosaic of randomly oriented modulated structures, very sim- 
ilar to that of the glass state. And moreover, if temperature 
is sufficiently small, this configuration develops an extended 
aging dynamics up to times far beyond any laboratory scale. 

We also found that, for any fixed magnetization, the h-T 
curve has a quadratic, growing functionality. In particular, 
very near the saturation, where the real system must pass from 
the modulated to the ferromagnetic state, the monotony of this 
functionality is a very important issue. It suggests a possi- 
ble interpretation based on the spin-bound for recently exper- 
imental observations on the phase diagram of ultrathin mag- 
netic films. lfl2ll Anyway, the fundamental impossibility of the 
classical spherical model to see the phase transition is also the 
responsible of the non-saturation of the system at any finite 
temperature. In this sense, a model carrying the nice prop- 
erties arising from the spherical restriction, but capable, as 
the Ginzburg-Landau, to account for the phase transition. ll22ll 
would be an important step forward. 

Such a model is introduced by means of the substitution 
of the global restriction by an infinite set of restrictions over 
finite size regions. The new dynamical equation is derived, 
and it is shown that this more complex but still linear equation 
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contains the classical spherical model as the limit case of para- 
magnetic high temperature state. It is suggested that, since the 
restriction is now over finite volumes in the modified model, 
once the correlation length becomes greater than this typical 
volumes size the system must fall into the Ising universality 
class. 
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Appendix A: Geometrical approach to some spatial correlations. 

Firstly we are intended to demonstrate, by means of geo- 
metric considerations, that spatial correlation C(R) in an en- 
semble of perfect striped configurations has an algebraic de- 
caying term of the form l/i?''*^^^/^, identical to that of ex- 
pression ( |39] l. One can perform such calculation by simply 
compute the spatial average 



C{R) = {(t>{x)4>{x + R)) 
with the order parameter taken as 

<f){x) = (j)o cos(fco • a; + $) 



(Al) 



(A2) 



d = 2 

Here we have to consider also the multiple angular orienta- 
tions of the modulated structures fco with respect to the axes 
origin. Now we have 



?) = 4>o cos(fco • 2^ + <&) 



(A6) 



Using polar coordinates we may write ko = {kg, a) and 
X = (r, 6). Now the random variables are the phase $ and 
the plane orientation of the structure a. In this way, the prob- 
ability of having a stripes configuration with values around 
certain $ and a is 



P($, a)d^ da 



da 
27r 27r 



(A7) 



And we have to integrate also the spatial average in polar co- 
ordinates. So, the calculation is 



C{R) 



leading to 



lim 



d$ r^"" da 

X 00 cos(fco ■ X + ^ + ka ■ R) 
X (j)Q cos(fco •£'+$) 



r dr d9 



r(R\ -^^^1 (h R\- '^0 cos(fcofl - 7r/4) 
C{R) - -MkoR) - Y 



(A8) 



(A9) 



where $ e [0, 27r] is included to represent the arbitrary rel- 
ative position between the modulated structure and the axes 
origin. 



d = 3 

Following the same analysis one starts now from an expres- 
sion of 6 



In one dimension the modulated structure can be written 



(j){x) — 00 cos(fcoa; + $) 



(A3) 



and the integration is simply over $ and x. Here, $ is the ran- 
dom variable containing all the relative displacements (phase) 
of the modulated structure with respect to the position x. The 
probability of finding a value between $ and $ + (i$ is equal 
to d(p/2TT. On the other hand, the spatial average is performed 
multiplying the spatial integral with limits ±L by a normal- 
ization factor 1/2L. So, we have to carry out the calculation 



2tt 



dx 



2L 



X 4>Q cos(fco2^ + $ + koR) 
X 00 cos(fcoa; + (f>) 



that leads to a perfect sinusoidal result 



C{R) = ^ cos(fcoi?) 



(A4) 



(A5) 



(j){x) — 00 cos(fco ■ X + ^) 



(AlO) 



in which ko and x are three-dimensional vectors thus rep- 
resenting alternating oriented laminae well-known in litera- 
ture. With a reasoning very close to that of d = 2, we use 
spherical coordinates and write vector as fco = (^Oi ot, ■0) ^nd 
X — (r, 6, If). The multiple orientations are given now by a 
and i/j and the stripes are characterized by fixing a and if) 
in such a way that 



P(<i>, a, il;)d^ da dip 



c?4> da sin(-0)d'0 
27r 27r 2 



(All) 



and taken the proper normalization for the spatial average in 
spherical coordinates 



CiR) 



lim 



2tt 



2^ 

L J2 



Jo 



sin(0) d'\\} da 

4 TT 



sin(6') dr dO dip 



10 Jo Jo 
X 00 cos(fco ■ x + ^ + ko ■ R) 
X 00 cos(fco • a; + $) 



(A12) 
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This leads to 



C{R) 



0Q sin(A:oi?) 0q cos(fcoi? — 7r/2) 



(koR) 



(A13) 



Finally, from expressions ( IA5I ). ( IA9I ) and dAOI l one can 
write, in the approximation kgR ^ 1, the following form of 
C{R) as function of the space dimension. 



C{R) 



t>l cos[fcofl- (ci~l)f] 
2 (fcoi?)^ 



(A14) 



and this is exactly the same functionality appearing in ( |39] | for 
spatial correlations, obtained in a completely different way. 



where it has been explicitly pointed the dependence of k^ and 
$ with X through squared brackets. This dependence is quite 
simple: while the modulus fco is constant over all the space, 
the relative phase $ and the rest of the d — 1 values Vli of 
the orientational components of k^ are constant only inside 
the limits of the finite regions and varies randomly between 
regions. So, we have to compute the average 

C{R) ^ 4>l (cos(A?o[x] • f + X (A16) 
X cos{kQ[x + R] ■ {x + R) + + R])) 

that is 



1. Glassy spatial correlation 

Lets analyze now the spatial correlation in a configuration 
of mosaics of randomly oriented modulated structures. In our 
approach we divide the space in squared d-dimensional vol- 
umes of equal linear size L inside which the local parameter 
is in a perfect striped pattern. The orientation of the modu- 
lated structures changes randomly from one area to the other. 
Under this condition one may write 

4i{x) = (/)o cos(A:o[^'] • X + (A15) 



C{R) (X J d'^x j d^'^-^'>n cos{kQ[x] ■ X + ^x]) x 

X cos{ko[x + R] ■ {x + R) + ^x + R]) (A17) 

One may note that, for R > L, that is, at distances higher 
that the typical domain size, the above expression becomes 
zero. This result is consistent with the exponential decay (l42l) 
in which the non-zero value of spatial correlations establishes 
up to certain length. For distances higher than this length the 
correlation suddenly vanishes. 



[1] P. M. Gleiser, F. A. Tamarit, S. A. Canas, and M. Montemurro, 

Phys. Rev. B 68, 134401 (2003). 
[2] E. A. Jagla, Phys. Rev. E 70, 046204 (2004). 
[3] L. Nicolao and D. A. Stariolo, Phys. Rev. B 76, 054453 (2007). 
[4] R. Di'az-Mendez and R. Mulet, Phys. Rev. B 81, 184420 (2010). 
[5] T. Garel and S. Doniach, Phys. Rev. B 26, 325 (1982). 
[6] C. Roland and R. C. Desai, Phys. Rev. B 42, 6658 (1990). 
[7] S. A. Pighm and S. A. Cannas, Phys. Rev. B 68, 134401 (2003). 
[8] R. Mulet and D. A. Stariolo, Phys. Rev. B 75, 064108 (2007). 
[9] R. E. Rosensweig, M. Zahn, and R. Shumovich, J. Magn. Magn. 

Mater. 39, 127 (1983). 
[10] M. Seul and R. Wolfe, Phys. Rev A 46, 7519 (1992). 
[1 1] M. Seul and D. Andelman, Science 267, 476 (1995). 
[12] N. Saratz, A. Lichtenberger, O. Portmann, U. Ramsperger, 

A. Vindigni, and D. Pescia, Phys. Rev Lett. 104, 077203 

(2010). 

[13] T. H. Beriin and M. Kac, Phys. Rev 86, 821 (1952). 

[14] R. J. Baxter, Exactly solved models in statistical mechanics. 

(Academic Press, London, 1982). 
[15] D. C. Brody, Phys. Rev E 49, 3665 (1994). 
[16] S. A. Cannas, D. A. Stariolo, and F. A. Tamarit, Physica A 294, 

362 (2000). 

[17] M. Grousson, G. Tarjus, and P Viot, Phys. Rev E 62, 7781 



(2000). 

[18] M. Grousson, V. Krakoviack, and G. Tarjus, Phys. Rev. E 66, 
26126 (2002). 

[19] R. Prozorov, A. F. Fidler, J. Hoberg, and P. C. Canfield, Nature 

Phys. 4, 327 (2008). 
[20] M. Tarzia and A. Coniglio, Phys. Rev Lett. 96, 075702 (2006). 
[21] R. Di'az-Mendez, A. Mendoza-Coto, R. Mulet, L. Nicolao, and 

D. A. Stariolo, arXiv:1011.5137vl, submitted to Phys. Rev B 

(2010). 

[22] A. Mendoza-Coto, R. Di'az-Mendez, R. Mulet, L. Nicolao, and 

D. A. Stariolo, to be published (2010). 
[23] L. F. Cugliandolo and D. S. Dean, J. Phys. A: Math. Gen. 28 

(1995). 

[24] M. Tarzia and A. Coniglio, Phys. Rev E 75, 01 1410 (2007). 
[25] C. Godreche and J. M. Luck, J. Phys. A: Math. Gen. 33, 9141 
(2000). 

[26] N. Saratz, A. Lichtenberger, O. Portmann, U. Ramsperger, 
A. Vindigni, and D. Pescia, Phys. Rev B 82, 184416 (2010). 

[27] P G. de Gennes, Solid State Commun. 10, 753 (1972). 

[28] B. Spivak and S. A. Kivelson, Phys. Rev B 70, 155114 (2004). 

[29] A. Abanov, V. Kalatsky, V. L. Pokrovsky, and W. M. Saslow, 
Phys. Rev B 51, 1023 (1995). 



